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Recently, M. Sieber and K. Richter achieved a break- 
through towards a proof of the BGS-conjecture by calculating 
a first semiclassical correction to the diagonal approximation 
of the orthogonal form factor for geodesic flow on a Riemann 
surface of constant negative curvature. 

In this note, we try to generalize the arguments. However, 
the solution proposed is not yet correct, because also other 
geometries must be taken into account. 



PACS numbers: 05.45.Mt, 03.65.Sq 



Introduction: Experimental and numerical work for many 
chaotic quantum systems has shown their spectral fluctu- 
ations to be faithful to random-matrix theory . How- 
ever, up to now, the equivalence of the spectral properties 
of random-matrix theory with the properties of classically 
chaotic quantum systems has only been a conjecture 
Following pioneering work of Sieber and Richter Q], we 
here give analytic evidence by a semiclassical treatment 
of geodesic flow on a surface of constant negative curva- 
ture with genus g > 2. The origin of universality of spec- 
tral fluctuations for chaotic quantum billiards and related 
systems with two-dimensional configuration space seems 
to lie in generic properties of self-intersections of periodic 
orbits in configuration space. Even though derived for a 
special system, our results could pave the way towards 
understanding universality: One would have to reveal the 
statistical properties of self-intersections as applicable to 
other dynamics. 

The so-called form factor, the Fourier transfom of the 
two-point correlator of the level density, is prediced by 
random-matrix theory for dynamics with time reversal 
invariance as JlJ 
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where r is the time in units of the Heisenberg time. In 
the interval < r < 1/2, one has the expansion K(t) = 

^2 ofc-fc 



2r 
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Let A 1 = T 7 /^/TrM 7 — 2 be the semiclassical ampli- 
tude of the classical periodic orbit 7 with period T 7 , ac- 
tion 5(7), and monodromy matrix M 7 . In the semiclas- 
sical framework, the form factor is then given by the fol- 
lowing double sum over periodic orbits Q 
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Here, Th denotes the Heisenberg time. According to 
Berry's 1985 insight [5j, one obtains the leading-order 
term 2r of random-matrix theory by retaining only the 
diagonal terms, 7 = 7', as well as the contructive in- 
terference of pairs of mutually time reversed orbits. We 
propose to establish the full r-expansion, of which the 
second order term — 2r 2 was found in the pioneering 
work of Sieber and Richter. They recognized that term 
as due to certain pairs of periodic orbits which are (in a 
certain sense exponentially) close neighbours in configu- 
ration space; one orbit in each pair has an additional self- 
intersection (with small angle e) relative to the other. A 
standard argument of Riemannian geometry ^ reveals 
the non-intersecting satellite orbit as unique. In sum- 
ming over the said pairs in equation ([!]), we shall be led 
to sums over intersections. 

Due to the mutual closeness of the orbits in each pair, 
their action difference S(j) — £(7') is not large compared 
to h; it is dominated by the behavior in the immediate 
neighbourhood of the intersection. In contrast to other 
" uncorrelated" pairs of orbits, the ones under consider- 
ation will therefore interfere constructively in the form 
factor. The pairs of orbits contributing the T 2 -term in 
K(t) involve one intersection. We shall show the r fc+1 
term to arise from families of orbit pairs involving k in- 
tersections (see Fig. 1). 

The model : As has been known for a long time, geodesic 
flow on a surface of constant negative curvature with 
genus g > 2 is chaotic 0. According to the Gauss- 
Bonnet theorem, a surface of constant negative curva- 
ture k = — 1 with genus g > 2 has area A = 47r(<? — 1). 
We are interested in self-intersections, because we sus- 
pect the universal behavior of classically chaotic systems 
to originate from the universal behavior related to the 
distribution of self-intersections. As has been known for 
a while, long closed orbits of length vT (only the product 
of velocity v and period T is of relevance) have a mean 
number of self-intersections fa] 
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We need to work with a distribution function P(e, T) 
for self-intersections with varying angle e (see Fig. 3), 
the epsilon integral of which is the above N(T), i.e. 
N(T) = 2 J +Tr P(e, T)de; here, a factor two arises since 
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Fig. 1 (a) 
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Fig. 2 (a) 
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Fig. 1 (b) 
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FIG. 1. (a) The "intersection box": When one orbit 
in a pair self-intersects, the companion avoids the crossing. 
There are two possibilities: Each arriving orbit may choose 
to self-intersect, (b) A k- family comprises 2 fc pairs of orbits 
with periods near T which have in the mean oc T 2 intersec- 
tions. The orbits in a pair differ only in the neighbourhood 
of precisely k points (the intersection boxes). 



the distribution function is even in e. Sieber and Richter 
showed that the order-T correction, not specified above, 
is responsible for the r 2 term in the form factor, while 
the order-T 2 term makes no contributions They 
have established the distribution in question for large T 
and in the interval < e < tt as 
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for negative e, P(e,T) is replaced by P{— e,T); c is a 
constant irrelevant for our discussion. 

We now turn to the action difference for a pair of or- 
bits arising from one intersection box. Following Richter 
and Sieber, we expand the action difference up to second 
order with respect to the intersection angle, 



AS(e,S 1 ,S 2 ) = ^{5 X + S 2 ) + 0(e 3 ) 
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Here, p — mv denotes the absolute value of the momen- 
tum, and Si + S 2 is the closest-approach distance of the 
non-intersecting orbit in the close encounter; that dis- 
tance is partitioned into 8\ and 82 by the point of in- 
tersection of the self-crossing partner orbit (see Fig. 3). 
Near the points of closest aproach, the outer orbit turns 
continuously from the stable direction of the left loop of 
the inner periodic orbit to the stable direction of the right 
loop. (Incidentally, this is an intuitive reason for which 
the outer orbit is periodic if the inner one is). 
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FIG. 2. (a) Allowed and forbidden encounters of orbits 
within a fc-family: A self-intersection is represented by a +, 
an avoided self-crossing by — . Allowed encounters within a 
fc-family are and (~) . Identical behavior of both orbits, 
(t) and (~), does not count as one of the k encounters defin- 
ing the fc-family. (b) Example for a sequence of possible en- 
counters in a 7-family. 
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FIG. 3. Geometry of the encounter of a self-intersecting 
orbits with its close-by partner 

In the system under discussion, for large times T 7 , the 
distances approach 8\ — 62 — * pe/ (2Xm) , where A is the 
Lyapunov exponent of the orbit. The Lyapunov exponent 
is the same for all orbits on a surface of constant negative 
curvature. In the limit of large T 7 , the action difference 
originating from one encounter becomes 
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that difference becomes small, eventually comparable to 
and even smaller than h, for e — > 0. And this is why pairs 
with small e make non-negligible "off-diagonal contribu- 
tions" to the form factor. At this point, we conclude our 
review of Refs. 

Multiple — loop contributions to the form factor : For 
the derivation of the small-r expansion, we will have to 
combine the contribution of all 2 fe different combinations 
of orbit pairs from the /c-families defined in Fig. 1 to get 
the term cx r k+1 . Fig. 1 shows two different periodic or- 
bits which are everywhere extremely close to one another 
except near k places in configuration space where one of 
them self-crosses, while the other one doesn't. In analogy 
to the somewhat related Hikami boxes p0[ of the theory 
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of disordered systems, we call these distinguished places 
in configuration space "intersection boxes". Let AS(ej) 
be the modulus of the action difference due to j-th inter- 
section box; it counts either with a positive or negative 
sign, to be denoted by sgii^ (7, 7') , depending on which 
of the two periodic orbits has the intersection. The to- 
tal action difference of the two orbits is proportional to 
the total length difference, which stems only from the 
intersection boxes, 



S(7)-S(7') = E s g n j(^') A5( £j ) 
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Now, consider a periodic orbit with period T. There 
are P(e\,T)dei possibilities for choosing an intersection 
with angle in the interval (ei,ei + de\). We assume dif- 
ferent intersection boxes to become uncorrelated for long 
periods T, because their number is much smaller than 
the total number of intersections of the periodic orbit in 
question. The joint distribution P(ei, £2, T) then 

simply involves the product of single-crossing distribu- 
tions P{cj,T). Since the intersections are ordered along 
the orbit, we must avoid overcounting by a factor 1/k in 
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The contribution to the term r k+1 is due to all 2 k 
orbit pairs of all £;-families. Let (—a, a) be the small in- 
tegration region for the intersection angles e, where the 
orbit pairs of interest exist. The semiclassical contri- 
bution of all orbit pairs from k intersection boxes in- 
volves a sum over all possible fc-families, and a sum 
over intra-family pairs, K^ h (r) oc £ k _ familios A 2 S(T - 

T l) Eintra-family pairs exp^Sfr) - tS(Y)). The Sum Over 

intra-family pairs requires integration over the k angles 
tj as well as summation, with the phase factor (6), over 
the 2 k possibilities of undergoing or avoiding crossings 
defined in Fig. 1(b), 

k — families 
p+a / 9 
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The first sum can be evaluated using the sum rule of 
Hannay and Ozorio de Almeida J|, £ 7 A 2 S(T-T 7 ) ~ T. 

Now, observe that in each factor lim^o J dejP{ej,T) 
exp(±ip 2 e 2 I (2mXK)) = U ± iN, there are two terms: A 
"non-universal" contribution N which is sensitive to the 



sign change in the exponent, and a "universal" contribu- 
tion U insensitive to that sign change. Since the orbit 
pairs exist only for small e, we have sine ~ e. After 
rescaling the integration variable as e = ep/ (y2mXh), 
and invoking the semiclassical limit of large %~ , the in- 
tegration range can be extended to (0,oo), whereupon 
we find 



U + iN^A dee (log e+...) e ±ir 
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The universal part U comes from the term —tt/2. The 
non-universal contribution has two components, one aris- 
ing from the logarithmic correction term in P(e, T) and 
another one, only represented by dots above, from the 
leading-order term oc T 2 in P(e,T). We thus get 



KU(t) oc E ( ] ) ( U + iN ) l ( U - iN ) k ~ l = ( 2U ) k ■ 
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Note that the non-universal part N (which, incidentally, 
depends on the Lyapunov exponent and would diverge 
semiclassically) cancels. Sticking in all factors, we finally 
obtain 
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It is reassuring to see a non-universal property like the 
Lyapunov exponent A cancel here. By summing over k 
we get 
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Of course, the foregoing semiclassical series sums up to 
the logarithm only in the range below half the Heisenberg 
time, i.e. for r < 1/2, where it gives the random-matrix 
result; for r > 1/2, the form factor is in principle deter- 
mined by its behavior up to r = 1/2 through unitarity 

By combining the present analysis with our previous 
work |nl], the symplectic form factor is derived as 
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Again, the series holds only up to half the Heisenberg 
time, which in the symplectic case means r < 1. Unitar- 
ity would have to be invoked to establish log(l — r) — > 
log 1 1 — \t\\, which random matrix theory claims for 
|r|<2. 

Conclusions: Obviously, the breakthrough achieved by 
Sieber and Richter has paved the way to a semiclassical 
understanding of the BGS-conjecture. But many open 
problems remain: From a physical point of view, the most 
intriguing one is the generalization beyond the geodesic 
flows considered here, in particular to ones with mixed 
phase spaces and to maps. Generalized time-reversal 
makes for another. Mathematicians will want to question 
much of the above, especially the assumed independence 
of intersections. 

Difficult as the open questions just mentioned may ap- 
pear, the very character of the above arguments suggests 
generalizability. First, we have already underscored the 
cancellation of the Lyapunov exponent from the univer- 
sal terms. No change is therefore expected for dynamics 
requiring a joint distribution P(e, A, T) and summation 
over A. Second, the existence of neighbouring pairs of 
unstable orbits for another type of billiards (hyperbolic 
boundary, flat space) was already confirmed by Sieber. 
The existence of a satellite to a self-intersecting orbit 
with small angle e follows from linerization of the flow 
around the latter unstable orbit and should hold even 
when stable orbits exist elsewhere. Third, maps in two 
dimensional phase spaces can be embedded in four di- 
mensional phase spaces as autonomous flows, at least if 
arising from Hamiltonians continuously and periodically 
modulated in time; projection on a two-dimensional sub- 
manifold should reveal self-crossing orbits and their satel- 
lites. Even the more familiar maps with Floquet oper- 
ators F — exp(—iH /h)exp(—iHi/h) should share the 
relevant properties since they can be thought of as gen- 
erated by Hamiltonians switching periodically between 
Hq and Hi; at the instants of switching, orbits in the 
extended phase space will still be continuous, but non- 
differentiable. The latter complication should not affect 
self-crossings in projections and the existence of non-self- 
crossing satellites. 

Support by the Sonderforschungsbereich "Unordnung 
und grofie Fluktuationen" of the Deutsche Forschungs- 
gemeinschaft is gratefully acknowledged. I enjoyed help- 
ful discussions with Martin Sieber, Uwe Abresch, Peter 
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Again, the series holds only up to half the Heisenberg 
time, which in the symplectic case means r < 1. Unitar- 
ity would have to be invoked to establish log(l — r) — > 
log 1 1 — |t||, which random matrix theory claims for 
|r|<2. 

Conclusions: Obviously, the breakthrough achieved by 
Sieber and Richter has paved the way to a semiclassical 
understanding of the BGS-conjecture. But many open 
problems remain: From a physical point of view, the most 
intriguing one is the generalization beyond the geodesic 
flows considered here, in particular to ones with mixed 
phase spaces and to maps. Generalized time-reversal 
makes for another. Mathematicians will want to question 
much of the above, especially the assumed independence 
of intersections. 

Difficult as the open questions just mentioned may ap- 
pear, the very character of the above arguments suggests 
generalizability. First, we have already underscored the 
cancellation of the Lyapunov exponent from the univer- 
sal terms. No change is therefore expected for dynamics 
requiring a joint distribution P(e, A,T) and summation 
over A. Second, the existence of neighbouring pairs of 
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The existence of a satellite to a self-intersecting orbit 
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them self-crosses, while the other one doesn't. In analogy 
to the somewhat related Hikami boxes [10] of the theory 
of disordered systems, we call these distinguished places 
in configuration space "intersection boxes". Let A5 , (e J ) 
be the modulus of the action difference due to j-th inter- 
section box; it counts either with a positive or negative 
sign, to be denoted by sgnj( 7 , 7 '), depending on which 
of the two periodic orbits has the intersection. The to- 
tal action difference of the two orbits is proportional to 
the total length difference, which stems only from the 
intersection boxes, 



5( 7 )-5( 7 ') = ^sgn j ( 7 , 7 ') ASiej) 



(6) 



Now, consider a periodic orbit with period T. There 
are P(ei,T)dei possibilities for choosing an intersection 
with angle in the interval (ei , e\ + de\ ) . We assume dif- 
ferent intersection boxes to become uncorrelated for long 
periods T, because their number is much smaller than 
the total number of intersections of the periodic orbit in 
question. The joint distribution P(ei, £2, e&, T) then 
simply involves the product of single-crossing distribu- 
tions P(ej,T). Since the intersections are ordered along 
the orbit, we must avoid overcounting by a factor 1/k in 



k 

P(e u e 2 ,...,e k ,T) = ^f[P(e j ) 



(7) 
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The contribution to the term r k+1 is due to all 2 k 
orbit pairs of all A;-families. Let (—a, a) be the small in- 
tegration region for the intersection angles e, where the 
orbit pairs of interest exist. The semiclassical contri- 
bution of all orbit pairs from k intersection boxes in- 
volves a sum over all possible fc-families, and a sum 
over intra-family pairs, K^ h (r) oc £ k _ families A 2 S(T - 
T i) Ei„tr a -f a mii y pairs exp(iS( 7 ) - iS( 7 '))- The sum over 
intra-family pairs requires integration over the k angles 
ej as well as summation, with the phase factor (6), over 
the 2 k possibilities of undergoing or avoiding crossings 
defined in Fig. 1(b), 
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The first sum can be evaluated using the sum rule of 
Hannay and Ozorio de Almeida [1], £ 7 A 2 S(T-T 7 ) ~ T. 

Now, observe that in each factor lim^o J " P( € j, T) 
exp(±ip 2 e 2 I (2mXh)) = U ± iN, there are two terms: A 



"non-universal" contribution N which is sensitive to the 
sign change in the exponent, and a "universal" contribu- 
tion U insensitive to that sign change. Since the orbit 
pairs exist only for small e, we have sine ~ e. After 
rescaling the integration variable as e = ep/ (\/2mXh), 
and invoking the semiclassical limit of large the in- 
tegration range can be extended to (0, oo), whereupon 
we find 



U + iNoc4 dee (log e> ...) e ±il 
Jo 

= ±ij o dp (log[pe ±i7r/2 ] + ...) e- p 

= -ir/2±i [J dp\og[p] e-' + ...) . (9) 

The universal part U comes from the term —ir/2. The 
non-universal contribution has two components, one aris- 
ing from the logarithmic correction term in P(e,T) and 
another one, only represented by dots above, from the 
leading-order term oc T 2 in P(e,T). We thus get 

KUW oc £ Q (U + iN)\U - iNf- 1 = (2U) k . 

(10) 

Note that the non-universal part N (which, incidentally, 
depends on the Lyapunov exponent and would diverge 
semiclassically) cancels. Sticking in all factors, we finally 
obtain 
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It is reassuring to see a non-universal property like the 
Lyapunov exponent A cancel here. By summing over k 
we get 

K 0Ith (T) =2t + tY, —j^ = 2r - r log(l + 2r). 

k=l 

Of course, the foregoing semiclassical series sums up to 
the logarithm only in the range below half the Heisenberg 
time, i.e. for r < 1/2, where it gives the random-matrix 
result; for r > 1/2, the form factor is in principle deter- 
mined by its behavior up to r = 1/2 through unitarity 

[!]■ 

By combining the present analysis with our previous 
work [11], the symplectic form factor is derived as 



OO 

^sympl.(r) = I + l£l- = I- I log(l - r). 

k=l 
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Fig. 1 (a) 



7 7' 



7 Y 



Fig. 2 (a) 

Y- 

7 - 



7 7' 



7 7' 



Fig. 1 (b) 
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FIG. 1. (a) The "intersection box": When one orbit in a 
pair self-intersects, the companion avoids the crossing. There 
are two possibilities: Each arriving orbit may choose to 
self-intersect, (b) A fc-family comprises 2 k pairs of orbits with 
periods near T which have in the mean oc T 2 intersections. 
The orbits in a pair differ only in the neighbourhood of pre- 
cisely k points (the intersection boxes). 



the epsilon integral of which is the above N(T), i.e. 
N(T) = 2 J v P(e,T)de; here, a factor two arises since 
the distribution function is even in e. Sieber and Richter 
showed that the order-T correction, not specified above, 
is responsible for the r 2 term in the form factor, while 
the order-T 2 term makes no contributions [4,9]. They 
have established the distribution in question for large T 
and in the interval < e < tt as 



P(e,T) 



TV 
2ttA 



4 log(e c) 
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for negative e, P(e,T) is replaced by P(—e,T); c is a 
constant irrelevant for our discussion. 

We now turn to the action difference for a pair of or- 
bits arising from one intersection box. Following Richter 
and Sieber, we expand the action difference up to second 
order with respect to the intersection angle, 



(4) 



Here, p = mv denotes the absolute value of the momen- 
tum, and Si + S 2 is the closest-approach distance of the 
non-intersecting orbit in the close encounter; that dis- 
tance is partitioned into Si and S 2 by the point of in- 
tersection of the self-crossing partner orbit (see Fig. 3). 
Near the points of closest aproach, the outer orbit turns 
continuously from the stable direction of the left loop of 
the inner periodic orbit to the stable direction of the right 



Fig. 2 (b) 
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FIG. 2. (a) Allowed and forbidden encounters of orbits 
within a fe-family: A self-intersection is represented by a +, 
an avoided self-crossing by — . Allowed encounters within a 
fc-family are and (~) . Identical behavior of both orbits, 
(+) and (~) , does not count as one of the k encounters defin- 
ing the /c-family. (b) Example for a sequence of possible en- 
counters in a 7-family. 
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FIG. 3. Geometry of the encounter of a self-intersecting 
orbits with its close-by partner 

loop. (Incidentally, this is an intuitive reason for which 
the outer orbit is periodic if the inner one is). 

In the system under discussion, for large times T 7 , the 
distances approach Si = S 2 — » pe/(2Xm), where A is the 
Lyapunov exponent of the orbit. The Lyapunov exponent 
is the same for all orbits on a surface of constant negative 
curvature. In the limit of large T 7 , the action difference 
originating from one encounter becomes 
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that difference becomes small, eventually comparable and 
even smaller than h, for e^O. And this is why pairs with 
small e make non-negligible "off-diagonal contributions" 
to the form factor. At this point, we conclude our review 
of Refs. [4,9]. 

Multiple — loop contributions to the form factor : For 
the derivation of the small-r expansion, we will have to 
combine the contribution of all 2 k different combinations 
of orbit pairs from the fc-families defined in Fig. 1 to get 
the term oc r k+1 . Fig. 1 shows two different periodic or- 
bits which are everywhere extremely close to one another 
except near k places in configuration space where one of 
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Recently, M. Sieber and K. Richter achieved a break- 
through towards a proof of the BGS-conjecture by calculating 
a first semiclassical correction to the diagonal approximation 
of the orthogonal form factor for geodesic flow on a Riemann 
surface of constant negative curvature. In the present paper 
we generalize and calculate both the orthogonal and symplec- 
tic form factors semiclassically. We find full agreement with 
random-matrix theory for times below half the Heisenberg 
time. 
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Introduction: Experimental and numerical work for many 
chaotic quantum systems has shown their spectral fluctu- 
ations to be faithful to random-matrix theory [1,2]. How- 
ever, up to now, the equivalence of the spectral properties 
of random-matrix theory with the properties of classically 
chaotic quantum systems has only been a conjecture [3]. 
Following pioneering work of Sieber and Richter [4], we 
here give analytic evidence by a semiclassical treatment 
of geodesic flow on a surface of constant negative curva- 
ture with genus g > 2. The origin of universality of spec- 
tral fluctuations for chaotic quantum billiards and related 
systems with two-dimensional configuration space seems 
to lie in generic properties of self-intersections of periodic 
orbits in configuration space. Even though derived for a 
special system, our results could pave the way towards 
understanding universality: One would have to reveal the 
statistical properties of self-intersections as applicable to 
other dynamics. 

The so-called form factor, the Fourier transfom of the 
two-point correlator of the level density, is prediced by 
random-matrix theory for dynamics with time reversal 
invariance as [1] 

f 2|r|-|r|In[l + 2|r|] : for |r| < 1 
^th(r)-| 2-Mlnf^i : for |r| > 1 , 

where r is the time in units of the Heisenberg time. In 
the interval < r < 1/2, one has the expansion K(t) = 

k=i 

Let Ay = Ty/^TrM-y — 2 be the semiclassical ampli- 
tude of the classical periodic orbit 7 with period T 7 , ac- 
tion 5(7), and monodromy matrix M 1 . In the semiclas- 
sical framework, the form factor is then given by the fol- 
lowing double sum over periodic orbits [1] 



^orth(r) = 

Um J- y a A*,S s ^ s ^l h 8{T - T t + Ty ) . (1) 
n^o T H ^ 7 7 v 2 ' v ' 

7,7' 

Here, Th denotes the Heisenberg time. According to 
Berry's 1985 insight [5], one obtains the leading-order 
term 2r of random-matrix theory by retaining only the 
diagonal terms, 7 = 7', as well as the contructive in- 
terference of pairs of mutually time reversed orbits. We 
propose to establish the full r-expansion, of which the 
second order term — 2r 2 was found in the pioneering 
work of Sieber and Richter. They recognized that term 
as due to certain pairs of periodic orbits which are (in a 
certain sense exponentially) close neighbours in configu- 
ration space; one orbit in each pair has an additional self- 
intersection (with small angle e) relative to the other. A 
standard argument of Riemannian geometry [6] reveals 
the non-intersecting satellite orbit as unique. In sum- 
ming over the said pairs in equation (1), we shall be led 
to sums over intersections. 

Due to the mutual closeness of the orbits in each pair, 
their action difference 5(7) —5(7') is not large compared 
to h; it is dominated by the behavior in the immediate 
neighbourhood of the intersection. In contrast to other 
"uncorrelated" pairs of orbits, the ones under consider- 
ation will therefore interfere constructively in the form 
factor. The pairs of orbits contributing the r 2 -term in 
K(t) involve one intersection. We shall show the T k+1 
term to arise from families of orbit pairs involving num- 
ber of k intersections (see Fig. 1). 

The model : As has been known for a long time, geodesic 
flow on a surface of constant negative curvature with 
genus g > 2 is chaotic [7]. According to the Gauss- 
Bonnet theorem, a surface of constant negative curva- 
ture k = — 1 with genus g > 2 has area A = Air{g — 1). 
We are interested in self-intersections, because we sus- 
pect the universal behavior of classically chaotic systems 
to originate from the universal behavior related to the 
distribution of self-intersections. As has been known for 
a while, long closed orbits of length vT (only the product 
of velocity v and period T is of relevance) have a mean 
number of self-intersections [8] 

2v 2 T 2 

N(T) = — r + 0(T) . (2) 

We need to work with a distribution function P(e,T) 
for self-intersections with varying angle e (see Fig. 3), 
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